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, We summarize recent ab initio studies of low-energy electroweak reactions of 

\ astrophysical interest, relevant for both big bang nucleosynthesis and solar neutrino 

■ production. The calculational methods include direct integration for up radiative 

and pp weak capture, correlated hyperspherical harmonics for reactions of ^ = 3, 4 
nuclei, and variational Monte Carlo for A = 6,7 nuclei. Realistic nucleon-nucleon 
and three- nucleon interactions and consistent current operators are used as input. 
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1 Realistic interactions and opportunities for astrophysics 

Combined partial- wave analyses [1] of essentially all elastic nucleon-nucleon 
(AW) scattering data have been used to construct highly accurate potential 
models [2,3,4] that describe the NN interaction with high precision, the "re- 
alistic" interactions. Available methods allow nuclear wave functions, binding 
energies, and electroweak matrix elements to be computed from these po- 
tentials for nuclei containing up to ten nucleons [5,6] with great success in 
matching the laboratory data. All of this is done with only bare interactions 
between pairs or triples of nucleons, and not with the effective interactions 
tailored to each class of problem that are more typical of nuclear physics. In 
particular, simpler potential models [7,8,9,10,11] of radiative captures in light 
nuclei need spectroscopic factors provided by experiment, while the resonat- 
ing group and similar techniques [8,12,13,14,15,16,17,18,19,20] rely on nuclear 
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interactions tailored to each system and use limited model spaces that affect 
the accuracy of cross section calculations in difficult-to-predict ways [19,20]. 

The mass and energy range accessible at present to calculations from re- 
alistic interactions and currents corresponds to two important areas of astro- 
physics: nucleosynthesis during the big bang [21] and hydrogen burning by 
the pp chain in low-mass stars including the sun [22,23]. The present paper 
reviews the methods and results of recent theoretical work on several rates 
important for these problems. 

The chemical composition of the universe just after the big bang was set 
by the freezeout of nuclear reactions when the universe was less than about 
five minutes old, a process called big-bang nucleosynthesis (BBN) [21,24]. The 
main products are hydrogen and ^He, in a ratio set mainly by weak interactions 
when the universe was about one second old. Small amounts of ^H, ^He and ''Li 
were also made, and their abundances after BBN provide crucial information 
on the mean baryon density in the universe [25]. Theoretical predictions of 
these abundances depend on the cross sections for eleven key nuclear reactions 
[26], all empirically determined quantities. At present, the reactions accessible 
to the methods discussed here are the radiative captures 7)^He, p{n,j)d, 
^He(Q;, 7)^Be, and ^H(q;, 7)'''Li. The first two of these are crucial for the 
yield of BBN, and the last three are crucial for the ^Li yield. 

The rates of energy and neutrino generation in the sun are set by the rate 
of the first reaction in the pp chain, the weak capture p + p — > + e"*" + z/e 
[22]. The cross section for this process is so small that it cannot be measured 
in the laboratory, so it must be provided by theory. The rates of other pp- 
chain reactions (reviewed by Haxton, Rolfs and Parker in this volume) do 
not appreciably affect energy generation in the Sun, but they do set the flux 
and energy distribution of solar neutrinos. In particular, the radiative captures 
^He(a, 7)'^Be and '^Be{p, 7)^B are accessible in the laboratory, but the Coulomb 
barrier and the low (~ 20 keV) energies involved make the cross sections very 
small and the experiments difficult. The absolute normalization is the crucial 
property of a cross section for astrophysics, and this is precisely what is most 
difficult to measure reliably in the laboratory. A theoretical treatment must 
be grounded in first principles in order to provide independent information 
about absolute cross sections. Theoretical calculations are also possible for 
weak decays of ^Be and ^B that actually produce the neutrinos. An additional 
process important for the solar neutrino spectrum, the weak capture ^He -|- 
p — ^ ^He + 6+ + (the "/iep" reaction), is accessible to these methods but 
not at all to experiment. 



2 The nucleeir Hamiltonians 

The Hamiltonians we use include nonrelativistic one-body kinetic energy, 
any of several modern, accurate two-nucleon potentials, and a three-nucleon 
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potential motivated by meson-exchange or chiral effective field theory: 



ijk ■ 



(1) 



i i<j i<j<k 



The kinetic energy operator is predominantly charge-independent (CI), but 
has a small charge-symmetry-breaking (CSB) component due to the differ- 
ence between proton and neutron masses. The modern high-precision two- 
nucleon (NN) potentials we consider are the Argonne vis (AV18) [2], CD- 
Bonn [3], Nijm-1, Nijm-Il, and Reid93 [4]. All these potentials are charge- 
dependent (CD), reproducing both the pp and np scattering lengths, as well 
as the deuteron binding energy, and fitting a total of 4301 NN scattering data 
from the 1993 Nijmegen partial wave analysis (PWA93) [1] up to Eiab — 350 
MeV with a /datum ^ 1. 

The AV18 potential is given by a sum of electromagnetic, onc-pion-cxchange 
(OPE), and short-range phenomenological terms. The strong interaction part 
is written as a sum of operator components: 



The first eight operators are the most important, reproducing the CI average of 
S- and P-wave phase shifts, while six additional terms quadratic in momentum 
are required to fit higher partial waves, and four small CD and CSB terms 
are required to explain fine differences between pp, np, and nn scattering. 
AV18 is local in every partial wave, and all partial waves are connected by 
this underlying operator structure. 

All the other potentials arc defined partial wave by partial wave. The 
Nijm-II and Reid93 potentials are local in each partial wave, while the Nijm- 
I central force includes momentum-dependent terms that give rise to non- 
local structures in configuration space. The CD-Bonn potentials have non- 
local components in both central and tensor parts of the interaction. All these 
modern potentials have about 40 parameters adjusted to fit the NN data. 
They cover a range of possible behaviors, as exhibited in their deuteron wave 
functions, shown in Fig. 1 of Ref . [27] . Despite these differences in detail, all 
five models give almost identical deuteron observables. 

The realistic NN potentials give varying amounts of binding in the triton, 
ranging from 7.6 McV for the local potentials to 8.0 MeV for the most non- 
local, compared to the experimental value of 8.48 MeV. In the ct-particle, 
the binding ranges from 24.0 to 26.3 MeV, compared to the experimental 
28.3 MeV [28]. Thus empirically, we need to add some three-nucleon (3A'") 
potential to obtain the proper binding and size of light nuclei. On the basis of 
meson-exchange theory we certainly expect 3A^ potentials to play a role, while 
chiral perturbation theory suggests there should be a rapid convergence in the 
contributions of many-body forces, so we may neglect four-nucleon forces. 
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The 3A^ forces we use are the Tucson-Melbourne (TM) potential [29], its 
chir ally- improved successor TM' [30], or one of the Urbana potentials [31]. 
These take the general form: 

Vijk = ^5'^ + V^[/ + V,% . (4) 

The TM and TM' potentials are constructed from two-pion exchange, where 
an intermediate tt/V scattering takes place in either S- or P-waves. The pion- 
range functions embedded in these terms have a cutoff that may be adjusted 
to fit the triton binding energy when used in conjunction with a particular 
AW potential. The Urbana potentials neglect the small S-wave term, but add 
a phenomenological short-range repulsion, V^^. In this case, the pion-range 
functions are taken to be the same as in the NN potential, while the strength 
of the terms are adjusted to reproduce the binding energy of the triton and 
to give a reasonable saturation density in nuclear matter. Urbana model IX 
(UIX) [32] was adjusted specifically to go with AV18. A more sophisticated 
set of Illinois model 3W potentials has been constructed recently [5], which 
include an additional three-pion ring term, VfffJ^^. The combined AV18/IL2 
Hamiltonian reproduces light nuclear binding energies very well up to ^4 = 10, 
but it has not yet been used in astrophysical rate calculations. 



3 Transition rates 

The capture processes of interest in the present study involve radiative or 
weak transitions between an initial two-cluster continuum state |p; JiMi, J2M2)^~^\ 
with clusters Ai and A2 having spins JiMi and J2M2, respectively, and relative 
momentum p, and a final 74-nucleon bound state \JaMa)- The cross section 
and polarization observables are obtained from matrix elements of the current 
operator connecting these two states. For example, the spin-averaged differ- 
ential cross section for radiative capture in the center-of-mass (CM) reference 
frame is written as 

\Ma M1M2 

where q is the momentum of the emitted photon, eA(q), A = ±1, are the 
spherical components of its polarization vector, and j(q) is the nuclear elec- 
tromagnetic current operator (see Sec. 4 below). The final bound state is 
recoiling with momentum — q. The p-direction has been taken to define the 
spin quantization axis, and p depends on the relative energy E=p'^ / {2/j), n 
being the A1-A2 reduced mass. The angle 6 is the angle between z and q, a 
is the fine-structure constant, uia is the bound-state mass, v = p/fi is the 
relative velocity, and the photon energy q is fixed by energy conservation. 
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The dependence of the obscrvables upon 6 can be derived from the expansion 
of the initial capture state |p; JiMi, J2^2)*'^'' into partial waves \LSJM)^~^\ 
and from the multipole expansion of the current operator [33] . The total cross 
section, however, is simply given by 



aiE) 



V l + g/mA(2Ji + l)(2J2 + l) 
X 5: ^ [ I Ei{q; LSJ) \^ + | Mi{q- LSJ) \^ ] , (6) 

1>1 LSJ 

where Ei and Mi denote the reduced matrix elements of the electric and mag- 
netic multipoles, i.e. Ei{q;LSJ) = (Ja || Ei{q) \\ LSJ)^~^^ and similarly for 
Mi{q;LSJ). 

The weak processes considered in the present study are proton weak cap- 
ture on and ^He, nuclear /3='=-decay and ^Be decay by electron (e-) capture. 
The cross section for the proton weak capture reactions is written as 

a^(E) = / 2.S(E. -E,)Ij: ^K/ I H„ | , (7) 



SeSv M1M2 



where H]^ is the weak-interaction Hamiltonian [34], E^ and Ef are the ini- 
tial and final state energies (including rest masses), and Pe and are the 
positron/electron and electron-neutrino momenta. The rates for the 0^ de- 
cays and e-capture are given by a similar expression as above, but with the 
flux factor 1/v removed and, in the case of the e-capture, the phase-space 
integration involves only the outgoing neutrino momenta. 

It is convenient to write the avereged weak-interaction Hamiltonian matrix 
element as E^es. E^vfo^K/ \Hw\i)? = {2nf L^r N"^ , where Gy is the 
Fermi coupling constant (G\/=1.14939x 10^^ GeV~^), and the lepton tensor 
L°"^ can be expressed in terms of the positron/electron and neutrino four veloc- 
ities [35]. The distortion of the charged- lepton wave function in the Coulomb 
nuclear field is also considered in L"'^ for processes with A > 4 [36]. The nu- 
clear tensor N'^^ can be expressed in terms of the reduced matrix elements of 
the Coulomb, longitudinal, electric and magnetic transition operators, corre- 
sponding to the multipole expansion of the nuclear weak current. Since the 
latter has scalar or vector (V) and pseudoscalar or axial-vector (A) compo- 
nents, each multipole consists of the sum of V and A terms, and the parity 
of Ith multipole ^-operators is opposite of that of Ith multipole A-operators. 
The parity of Coulomb, longitudinal, and electric Ith multipole V-operators 
is (— )^, while that of magnetic Ith multipole ^-operators is (— )^"'"^. 

Finally, in later sections, the results for capture reactions between charged 
clusters will be presented in terms of the astrophysical S'-factor, related to the 
total cross section <j{E), via S{E) = Ea{E)exp{2'n' ZiZ2a/v), where Zi and 
Z2 are the atomic numbers of the incoming clusters. The exponential factor is 
the Coulomb penetration factor. 
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4 The nucleeir transition current operators 

The model for the nuclear electroweak current has been most recently 
reviewed in Refs. [33,35]. The current is written as sum of one- and many- 
body components that operate on the nucleon degrees of freedom. The one- 
body operators are derived from the nonrelativistic reduction of the covariant 
single-nucleon electroweak currents and include terms up to order where 
m is the nucleon mass [33,35] . Inclusion of only one-body terms is known as the 
"impulse approximation" (lA). The two-body operators are discussed briefly 
in the following two subsections. 

4-1 The electromagnetic two-body currents 

The two-body current has "model-independent" and "model-dependent" 
components, in the classification scheme of Riska [37]. The model-independent 
terms arc obtained from the two-nucleon interaction and arc constructed [38] 
so as to satisfy current conservation. The leading operator is the isovector "tt- 
like" current obtained from the isospin- dependent spin-spin and tensor interac- 
tions [33] . The same components of the interaction also generate an isovector 
"p-like" current. Additional model-independent isoscalar and isovector cur- 
rents arise from the isospin-independent and isospin-dependent momentum- 
dependent interactions. However, these currents are short-ranged and numer- 
ically far less important than the vr-like current. 

The model-dependent currents are purely transverse and therefore cannot 
be directly linked to the two-nucleon interaction. Present results include the 
contributions due to the isoscalar pnj and isovector cuttj transition currents 
as well as those due to the isovector current associated with excitation of 
intermediate A resonances. Two different approximations have been used to 
treat the A degrees of freedom [33]; one is based on first-order perturbation 
theory, using the static A approximation, while the other is the transition- 
correlation-operator (TCO) scheme [39], essentially a scaled-down version of 
the full A^+A coupled-channel method. Comparisons between results obtained 
within these two approximate schemes have been reported in a number of 
studies [39,40,41]. 

The presence of 3A" interactions in the Hamiltonian requires corresponding 
three-body currents in order to fulfill current conservation, but these were 
found to be numerically unimportant in studies of the trinucleon magnetic 
form factors at intermediate momentum transfer [41]. 

4-2 The weak two-body currents 

The nuclear weak current consists of vector {V) and axial- vector (A) parts, 
with corresponding one- and two-body components. The weak vector current 
is constructed from the isovector part of the electromagnetic current, in ac- 
cordance with the conserved- vector-current hypothesis [34]. 
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Some of the two-body axial-current operators are derived from tt- and 
p-meson exchanges and the pvr-transition mechanism. These mesonic opera- 
tors have been found to give rather smaU contributions in weak transitions 
involving few-nucleon systems [27,35,42]. The two-body weak axial-charge op- 
erator includes a pion-range term from soft-pion theorem and current algebra 
arguments [43], as well as short-range terms associated with scalar- and vector- 
meson exchanges. The latter are obtained consistently with the two-nucleon 
interaction model, following a procedure [44] similar to that used to derive the 
corresponding weak vector-current operators [35]. 

The dominant two-body axial current operator, however, is that due to 
A excitation [27,35]. Since the A^A transition axial- vector coupling constant 
g\ is not known experimentally, the associated contribution suffers from a 
large model dependence. To reduce it, g\ has been adjusted to reproduce the 
experimental value of the Gamow- Teller (GT) matrix element in (3 decay 
(cf. Table 6). As in the case of the electromagnetic current, the perturbative 
and TOO approximations have been used to calculate the A-excitation ax- 
ial current contributions. We emphasize that results obtained within the two 
schemes are typically very close to each other once g\ is fixed, for each scheme, 
to reproduce the GT matrix element (cf. Table XV of Ref. [35]). 

Some of our weak transition studies in few-nucleon systems have also been 
carried out by using weak current operators derived within an effective-field- 
theory (EFT) approach, in which pions and nucleons are the explicit degrees of 
freedom and terms up to next-to-next-to-next-to-leading order (N'^LO) are in- 
cluded. Explicit expressions for these EFT operators can be found in Ref. [45]. 
Here we only note that: i) the two-body operators are regularized by a cut- 
off that defines the energy /momentum scale of EFT below which the chosen 
explicit degrees of freedom are valid; ii) the axial-current terms depend on an 
unknown parameter d^, that gives the strength of a counterterm and is fixed, 
for a given value of the cutoff, by again reproducing the '^H GT matrix ele- 
ment. The premise of EFT is that physical observables should not depend on 
the cutoff as long as the latter is chosen in a physically reasonable range. Such 
an expectation is indeed borne out in the EFT calculation of the pp capture 
discussed below. 

5 Wave functions 

Wave functions (bound and continuum) for the present calculations are 
computed by two means: correlated hyperspherical harmonics (CHH) and 
variational Monte Carlo (VMC). In the CHH method, the nuclear wave func- 
tion is expanded on a suitable basis and the unknown expansion parame- 
ters are calculated using appropriate variational principles. This technique is 
very successful in describing nuclear ground states with A < 4 [46,47,48,49], 
the Nd scattering states both below and above the deuteron breakup thresh- 
old [50,51,52,53,54], and the p^He and n^B. systems below the ^He and ^H 
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breakup thresholds [55,56]. VMC is an approximate method that uses Monte 
Carlo sampling to perform numerical quadratures and obtain upper bounds 
to ground and excited state energies for 3 < ^4 < 10 nuclei [6,57]. There is also 
the Green's function Monte Carlo (GFMC) method that starts from a VMC 
trial function and employs Monte Carlo integration to evaluate an imaginary- 
time path integral that projects out exact states of a nucleus. Extension of 
the Monte Carlo methods to describe continuum states is more difficult, as 
discussed below, and the GFMC method has not yet been used to compute 
reaction cross sections. 

5. 1 Hyperspherical harmonic wave functions for bound states 

We first consider the CHH method applied to the A = 3 bound state. The 
wave function ^ of a three-nucleon system with total angular momentum JM 
and total isospin TT^ can be written as [46] 

* = E«a-n(P)E^a,n(p) • (8) 
a,n p 

Here p is the hyper-radius, p = {x"^ + yf)^^"^, with the Jacobi coordinates Xj 
and Yi defined respectively as Xj = — and yj = (r^ -|- — 2 Fj) / ^/S. The 
index p denotes an even permutation of the particle indices i,j,k. The functions 
Ha,n are antisymmetric under the exchange j ^ k. They are given by the 
product of an angle-spin-isospin function and a hyperspherical polynomial 
Pa,n{4'i)i (pi being the hyper-angle defined as (pi = cos^^{xi/ p). The index a 
specifies the set of the spectator i and pair jk orbital angular momenta, la 
and La, spin, and isospin, which are coupled to produce the given JMTT^ 
quantum numbers and the appropriate parity. The Pa,n{4'i) is proportional to 
a Jacobi polynomial; the second index n runs over all non-negative integers 
and specifies the order K = + La + 2n of the polynomial. The product 
of the Jacobi polynomial with the spherical harmonics Yi^{Sti) and yL„(yi) 
contained in is by definition a hyperspherical harmonic function. Finally, 
correlation factors, which arc functions of the jk relative distance, are included 
in Ha,n to account for the strong state-dependent correlations induced by 
the NN interaction. These correlation functions are solutions of suitable two- 
body Schrodinger equations [46] . Their presence greatly improves the rate of 
convergence of the expansion of Eq. (8); the results quoted in this work have 
been obtained with 23 angle-spin-isospin channels, for which the maximum K 
is 16 [48]. To obtain a comparable accuracy without correlation factors, the 
maximum value of K needs to be about 180 [58]. 

The Rayleigh-Ritz variational principle is used to determine the hyper- 
radial functions Ua,n{p) s-nd the ground state energy E. Carrying out the 
variation with respect to the functions Ua,n{p), and performing the spin- 
isospin sums and the integration over the angular and hyper-angular variables 
yields a set of coupled second-order differential equations for the Ua,n{p)- This 
is then solved numerically [46,50]. 
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Table 1 

Calculated AV18/UIX and experimental binding energies of light nuclei in MeV. 



'^Z{J'';T) CHH VMC(I) VMC(II) GFMC Expt 





8.479 


8.227(6) 




8.461(6) 


8.482 


'He(l+;i) 


7.750 


7.476(6) 




7.708(6) 


7.718 


%e(0+;0) 


27.89 


27.40(3) 




28.31(2) 


28.30 


6Li(l+;0) 




28.05(5) 


28.16(5) 


31.25(8) 


31.99 


^Li(i-;i) 




33.07(7) 


32.47(7) 


37.5(1) 


39.24 


7Li(i-;i) 




33.13(7) 


32.23(7) 


37.6(1) 


38.76 






31.49(7) 


30.60(7) 


35.9(1) 


37.60 


'Be(r;^) 




31.58(8) 


30.71(7) 


35.9(2) 


37.17 



The extension of the CHH approach to the study of the a-particle is con- 
ceptually straightforward but numerically much more involved [47,55] . In anal- 
ogy with the three-body case, the wave function is written as in Eq. (8) and 
the Rayleigh-Ritz variational principle is again used to determine the hyper- 
radial functions and ground state energy E. However, it should be noted that 
in the case of A = 4, there are two sets of Jacobi coordinates, correspond- 
ing to the 1-1-3 and 2-1-2 partitions, and both partitions have been used in 
the definition of the angle-spin- isospin functions Y^. Furthermore, there are 
three hyperspherical coordinates, the hypcr-radius and two hyper-angles, and 
only one hyper-angle depends on the partition considered. The hyperspheri- 
cal polynomial of Eq. (8) depends in this case on two non-negative integers, 
which specify the order of two individual Jacobi polynomials, functions of the 
hyper-angles. Finally, the correlation factors are products of pair-correlation 
functions and are obtained by the same procedure as in the three-body case. 
More details are given in Ref. [47]. 

The results for the binding energies of ^H, ^He and "^He obtained with the 
CHH method for the AV18/UIX Hamiltonian are given in Table 1 [35,48]. Note 
that an accuracy of 1 (10) keV can be achieved for the three-body (four-body) 
problem. 

5.2 Variational Monte Carlo wave functions for bound states 

A good VMC wave function, ^y( J'^; A, T), is constructed from products of 
two- and three-body correlation operators acting on an antisymmetric Jastrow 
wave function with the appropriate quantum numbers: 

. (9) 



I* 



V 



TNI 
ijk 



i<j<k<A 



3 n ( 

i<j<A 



1 + Uij] 



9 



Here S denotes symmetrization over the noncommuting pair correlation opera- 
tors Uij, which inchidc spin, isospin, and tensor terms, Uij = J2p=2,6'^p{f^ij)Oij, 
and Of, are the same static operators that appear in the AW potential, Eq.(2). 
The U^jk^ are correlations induced by the three-nucleon interaction, and have 
the same spin-isospin structure as the Vijk of Eq.(4). 

The Jastrow wave function, ^j, for s-shell nuclei has the simple form 

^a{JMTT,)) . (10) 

Here fci^ij) and ff^f^ are central two- and three-body correlation functions 
and $^ is a Slater determinant in spin-isospin space, e.g., for the a-particle, 
|<l>4(0000)) = A\p ] p i n ^ n I). The functions /c(?") and Up{r) are gener- 
ated by solving a set of coupled differential equations containing the bare NN 
potential with appropriate asymptotic boundary conditions [59]. Variational 
parameters in these equations arc adjusted to minimize the energy expecta- 
tion value, Ev = {^v\H\^v), which is evaluated by Metropolis Monte Carlo 
integration [60]. 

The ^ J for p-shell nuclei is more complicated, with single-nucleon p-shell 
orbitals and LS couphng to obtain the desired JM value [61]. We sum over 
all allowed spatial symmetries [n] [62] and allow for three types of central pair 
correlations depending on which shell the nucleons are in: 

j<j<fe<4 i<j<4 LS[n] k<A<l<A 

n /pp^'"V/m)I^A(L5[n]JMTT,)i234:5...A)) } ■ (H) 

i<l<m<A 

The operator A indicates an antisymmetric sum over all possible partitions 
into 4 s-shell and [A — 4) p-shell particles. The pair correlation for both par- 
ticles within the s-shell, fss-, is similar to the of the a-particle. The pair 
correlations for both particles in the p-shell, f^^"'\ and for mixed pairs, fgp^^"'\ 
are similar to fss at short distance, but their long-range structure is adjusted 
to give appropriate clustering behavior, and they may vary with LS[n\. The 
single-particle wave function components are: 

\^A{LS[n]JMTT,)u34:5...A) = | $4 (0000) 1234) x | [] 'Z'^''^"' (^aO (12) 

4<1<A 

\ [ n yimi{^al)]LMj^ln] X [ [] XK^^s)]sMs \ X [ [] i^/(^^t)]rT,) ■ 
[ 4<1<A 4<1<A ^ ) jj^ 4<1<A ^ 

The (pp^^^^ (Rai) are p-wave solutions of a particle in an effective a-N potential 
that has Woods-Saxon and Coulomb parts. They are functions of the distance 
between the center of mass of the a-core and nucleon /, and may vary with 
LS[n]. The wave function is translationally invariant, so there is no spurious 





n fijk 




n fcir^,) 




i<j<k<A 




i<j<A 
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center of mass motion. 

Two different types of have been constructed in recent VMC calcula- 
tions of light p-shell nuclei: Type I, which is a shell- model trial function [63], 
and Type II, which is a cluster-cluster trial function [64,65]. In Type I trial 
functions, the ^^'^["'(r) has an exponential decay at long range, with the depth, 
range, and surface thickness of the Woods-Saxon potential serving as varia- 
tional parameters. The /^^'^'"^ go to a constant near unity at long range, while 
the fpp^^^ have a small long-range tail that is larger for states of lesser spatial 
symmetry [n] [63]. In Type II trial functions, ^^'^["■^(r) is again the solution of 
a p-wave Schrodinger equation with a potential containing Woods-Saxon and 
Coulomb terms, but with an added Lagrange multiplier that turns on at long 
range and imposes the boundary condition [^^'^["'(r — > oo)]" oc W^fc^(27r)/r. 
Here Wkmi^^f ) is the Whittaker function for bound-state wave functions in a 
Coulomb potential and n is the number of p-shell nucleons. The cluster sep- 
aration wave number 7 is taken from experiment. The accompanying f^p^""^ 
goes to unity while the f^^'^^ are taken from the exact deuteron wave func- 
tion in the case of ^Li, or the VMC '^H ('^He) in the case of ^Li (^Be). 
Consequently, Type II trial functions factorize at large cluster separations as 
iJJa'^JrWkmi'^iraT) /rar, whcrc tpa and ipr are the wave functions of the 
clusters and far is the separation between them [64,65]. 

For either type of trial function, a diagonalization is carried out in the 
one-body basis to find the optimal values of the PlsIu] niixing parameters for 
a given [J'^; T) state. Current best VMC energies for nuclear states of astro- 
physical interest are given in Table 1 for the AV18/UIX Hamiltonian along 
with CHH and GFMC results and experimental values. These VMC results are 
about 1 MeV more bound for ^4=6, 7 nuclei than previously reported because 
of recent improvements in the parametrization of the trial functions; however, 
results for various transition matrix elements are not significantly altered. The 
VMC energies for A=3,4 nuclei are 1-2% above CHH or GFMC results, but 
are 10-15% higher for A=6,7 than GFMC results, which are themselves 5% 
high compared to experiment with this Hamiltonian. 



5.3 Hyperspherical harmonic wave functions for continuum states 

The CHH method has also been applied to the study of scattering prob- 
lems. In particular, pd and nd systems have been studied both below [50,51] 
and above [52,53,54] the deuteron breakup threshold. The p^'Ylc and n^H sys- 
tems have also been investigated below the "^He (^H) breakup threshold [55,56]. 
We now review how the CHH method works for the N + A scattering below 
the mass-A nuclear breakup threshold. 

The wave function ^^r+'^, having incoming orbital angular momentum L 
and spin S coupled to total JM, is written as the sum of a core function \1';^^ 
and an asymptotic function ^^•S'-^^. xhe function describes the system 
in the region where the particles are close to each other and their mutual 
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Table 2 

Calculated A = 3 doublet (a^) and quartet (a^), and A = A singlet (og) and triplet 
(at) scattering lengths in fm [40,55]. Experimental results are from Refs. [67,68]. 





nd 


pd n^H p^Re 




o-d aq 


ad aq as at ag at 


CHH 


0.63 6.33 


-0.02 13.7 4.32 3.80 11.5 9.13 


Expt 


0.65(4) 6.35(2) 


10.8(2.6) 8.1(5) 



interactions are large, and vanishes in the limit of large inter-cluster separation. 
It is obtained by an expansion on the same CHH basis as discussed above for 
the bound state. In the asymptotic region the inter-cluster interactions are 
negligible, and •^^^■^^ is written as (in the p + A case) 



^a^-""" = E E h ® Ms' ® l^L'(fM) 



i L'S' 
X 



JM 



OlL'OsS' 7- ^ ^LS,L'S'W — 7- 9{rpA) 

HI pA 1^1 pA 



(13) 



Here (pA, fpA and k are the mass-74 cluster wave function, the proton and 
74-cluster relative distance, and the magnitude of the relative momentum, 
respectively. The functions and Gl are the regular and irregular Coulomb 
functions, respectively. For n+A scattering, Fl{x)/x and Gl{x)/x are replaced 
by the regular and irregular spherical Bessel functions. The function g{rpA) 
regularizes G^^krpA) at rpA — 0, and g{TpA) — 1 as rpA > 10 - 12 fm. Finally 
Rls L'S'i^) are the i?-matrix elements that determine phase shifts and (for 
coupled channels) mixing angles at the energy k'^/{2fi), fi being the N + A 
reduced mass. The sum over L'S' includes all values compatible with a given 
total J and parity. Both the matrix elements R'ls,L'S'i^) hyper-radial 
functions occurring in the expansion of are determined by applying the 
Kohn variational principle. 

The nd, pd, and p^He scattering lengths predicted by the AV18/UIX model 
are shown in Table 2 [40,55]. Also listed in the table are the n^H scattering 
lengths predicted by the older AV14/UVIII Hamiltonian [59,66]. There is ex- 
cellent agreement between the theoretical predictions and the available data. 



5.4 Phenomenological continuum states for A> 4 



Few calculations of continuum states consisting of two nuclei have been 
done in the quantum Monte Carlo formalism, and several difficulties in these 
calculations remain to be addressed. We expect these problems to be tractable, 
and work on them continues. In the meantime, we have adopted a semi- 
phenomenological prescription for the continuum-state wave functions that 
enter into capture calculations for the A > 4 systems [64,65]. 

At large separations between nuclei in a continuum state, the wave function 
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is the product of the wave functions 0i and ^2 of the two nuclei with a two- 
body correlation that describes the (nuclear plus Coulomb) scattering of the 
two nuclei from each other. This is 



^LSJM 



[01«)02]5®>L(ri2) 



J JM 



— h tan dLSj[k) — ; 



,(14) 



where ri2 is the separation between the initial nuclei. This is analogous to 
Eq.(13) without coupled channels, and the angular momentum labels have 
the same meanings. 

To obtain approximate continuum wave functions, we have adopted the 
form 



JM 



for the wave function. The operator A antisymmetrizes the wave function un- 
der exchange of nucleons between the two (internally antisymmetric) clusters. 
The short-range Gij operator correlations between nucleon pairs in different 
clusters produces distortions in the nuclei when they occupy the same space, 
and they go to unity at less than 2 fm interparticle separation. We used corre- 
lations Gij found in nuclear-matter calculations [69] and found that they had 
essentially no effect on the cross sections. 

We generated the radial function ■0i2^'^('^i2) from a Schrodinger equation 



2/1 
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+ 145j(ri2) \ Mil') = E{T^2i>\l'). (16) 



The appropriate reduced mass for the nuclei is ^lyi- The effective potential 
function VL5j(r 12) was chosen to be the sum of a short-range nuclear inter- 
action (Woods- Saxon or similar, with central and spin-orbit terms) and the 
Coulomb interaction. We chose Vlsj{jv2) to reproduce the experimentally de- 
termined phase shifts SLSj{k)- For both A = 6 and A = 7 a-captures, several 
suitable choices of Vlsj were already available in the literature, so we chose 
from among those [7,8,9,10,70,71,72], making minor adjustments to reproduce 
the well-known width and location of the 3^ resonance in ®Li. The dependence 
of our results on the different Vlsj indicates the amount of uncertainty intro- 
duced by this approach. To satisfy the Pauli principle, the functions Vlsj need 
deep potential wells with forbidden states to produce the correct number of 
nodes in ipi-f'-' [73], as well as parity dependence. The ^Li Vlsj also contains 
a spin-orbit term to describe three low-lying resonances important for the ra- 
diative capture. It can in principle also admit a tensor interaction (requiring 
modification of Eqs.(14) and (16)), but it is not unambiguously required by 
the data, and we did not include it. 

The J'^ = 1^ wave functions in ad should also be to be orthogonal to 
the ^Li ground state, but there is no guarantee of this in the present approach 
because "012^"^ was not obtained from the underlying A^'A'" interaction. We found 
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that this property is only important for magnetic dipole transitions proceeding 
from S waves. For this matrix element, we kept only the part of "if^f^l^^ 
orthogonal to the VMC ground state. 

Finally, the states 0i and 02 appearing in Eq.(15) were computed directly 
from the two- and three-nucleon forces by the VMC method, except that the 
deuteron wave function 0^ was a direct numerical solution of the two-body 
Schrodinger equation for the deuteron. 

6 Results 

6.1 The p(n,j)d radiative capture below 1 MeV 

Historically, the radiative capture of thermal neutrons on protons played 
a crucial role in establishing the quantitative importance of two-body current 
effects in photonuclcar obscrvables [74]. Their inclusion resolved the long- 
standing discrepancy between the calculated impulse approximation (lA) cross 
section and the measured value. In this section, we briefly summarize the 
results of cross section calculations for this process from thermal energies up 
to energies of ~ 1 MeV; the energy range between a few tens of keV up to 200 
keV is particularly relevant for primordial nucleosynthesis. 

At thermal energy, the cross section a'^ is due to the well-known isovec- 
tor magnetic dipole (Mi) transition connecting the ^Sq np scattering state 
to the deuteron bound state. There is in principle an isoscalar Mi transition 
proceeding through the ^Si np continuum state, but the associated contri- 
bution (at thermal energy) is strongly suppressed in lA by orthogonality; 
furthermore, isoscalar two-body currents, such as those originating from the 
momentum-dependent components of the two-nuclcon interaction or from the 
pK^ transition mechanism, are numerically far less important than isovector 
currents, and play a very marginal role in this transition. 

The calculated values for a"^ are listed in Table 3, both for one-body cur- 
rents alone and for the one- and two-body currents. The lA results are to a 
large extent determined by the ^5*0 scattering length and deuteron binding 
energy, which the above interactions reproduce exactly. Indeed, in effective 
range theory, one finds that a'^ ~ 300 mb. Among the two-body current con- 
tributions, which collectively amount to approximately 9% of the total cross 
section, the largest one, about 60%, comes from the currents associated with 
pion exchange. The next to largest, about 25%, is from currents involving A 
excitation, that have been treated in perturbation theory here. 

As the energy increases from a few to several hundreds of keV, contribu- 
tions to a'^ from the electric dipole (Ei) transitions proceeding through the 
^Pj np scattering states become more important, and in fact dominant at 
energies larger than 300 keV. The Siegert form is used for the Ei operator, 
and relativistic corrections to the charge density, such as those associated with 
spin-orbit and pion-exchange contributions [76,64], have been neglected. 
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Table 3 

Total cross-section a"^ in mb for thermal np radiative capture; data from Ref. [75]. 



NN model 


AV18 


Nijm-I 


CD-Bonn 


Expt 


One-body 


304.6 


305.4 


306.5 




Full 


334.2 


332.5 


331.6 


332.6(7) 



The results of calculations based on the AV18 and CD-Bonn interactions 
are shown in Fig. 1. As the energy increases the model dependence, already less 
than 1% at thermal energies, becomes negligible at 1 MeV, where the Ei con- 
tribution dominates. This is to be expected, since the Ei transition is mostly 
sensitive to the long-range part of the wave functions. Fig. 1 demonstrates that 
the model dependence of the predictions in the conventional approach based 
on realistic interactions is considerably less than 1% in the energy range rele- 
vant for primordial nucleosynthesis. It could be further reduced by adjusting, 
for example, the transition magnetic moment /i-yATA of the A-excitation cur- 
rent to precisely fit the thermal neutron capture cross section. In all results 
reported so far /i^ata has been taken to be 3 nm, a value consistent with an 
analysis of data at resonance. 

Finally, it is worth emphasizing that the model dependence referred to 
above is comparable to the theoretical uncertainty inherent in calculations 
based on effective field theory [77,78]. These calculations involve contact op- 
erators whose strength has to be constrained by experimental data. In the 
more recent study of Ref. [78], for example, there is a contact four-nucleon- 
photon operator, whose Mi (Ei) strength is fixed by the measured thermal 
capture cross section (deuteron photo-disintegration data at threshold). 

8 
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Fig. 1. Predictions for np radiative capture cross section, calculated with the AV18 
and CD-Bonn interactions. Also shown for AV18 is the Mi contribution to a"' . 
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Table 4 

Square of the overlap integral A{E = 0) for five modern AW interaction models. 



NN model 


AV18 


Nijm-I 


Nijm-II 


Reid93 


CD-Bonn 


(one-body) 


6.965 


6.965 


6.971 


6.974 


6.985 


A2 (full) 


7.076 








7.060 



6.2 The pp weak capture 



The theoretical description of the proton weak capture on protons was 
first given by Bethe and Critchfield [22], who showed that the associated rate 
was large enough to account for the energy released by the sun. Since then, a 
series of calculations have refined their original estimate either by computing 
the required wave functions more accurately [79,80,81,82] or by using more 
realistic models for the nuclear transition operator [83.84.85]. However, the 
most recent studies [27,45] offer an integrated treatment of these two aspects 
with emphasis on a reliable estimate of their associated uncertainties. 

The pp weak capture is induced by the weak axial current operator, and 
the ^'-factor is conventionally written in terms of a dimensionless parameter 
A{E) [27], where E is the pp relative energy. When axial two-body currents 
are neglected, it is given by 

oo 

A{E) = (7V2)'^' ^TT^ / dru{r)xo{r;k) , (17) 

where 7 is the deuteron wave number, k is the pp relative momentum [E=k'^ /nip), 
and Co{k) is the Gamow penetration factor. In the overlap integral, u{r) 
is the S-wave component of the deuteron wave function and Xo(^:^") is the 
pp ^Sq continuum wave function. The ^5*0 wave function is the solution of a 
Schrodinger equation with a pp interaction including, in addition to the nu- 
clear and Coulomb terms, terms of second order in the fine structure constant. 
These are the vacuum polarization and two-photon exchange terms, as well 
as corrections originating from the finite size of the proton [27]. Properly ac- 
counting for these additional electromagnetic interaction components leads to 
some practical problems in integrating the Schrodinger equation, matching 
the Coulomb functions, and extracting the phase shift. The most important 
correction to the Coulomb interaction between protons is the vacuum polar- 
ization, which reduces the cross section (proportional to A^) by about 1%. 
Other fine details of the electromagnetic interaction increase the cross section 
by about 0.1%. This is in part canceled by a net 0.03% reduction in cross 
section from the correct relativistic treatment of the deuteron wave number, 
7. Including just the axial one-body operator, five modern interaction models 
differ by only 0.3% in the calculated cross section as illustrated in Table 4. 

The model for the axial two-body current contains the simplest possible 
two-body operators that give an adequate description of the longest-range 
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mechanisms. The overall strength of the leading operator due to A excitation 
is adjusted as described in Sec. 4.2. The contributions due to exchanges of 
heavier mesons, such as the Ai [86], or renormalization effects arising from 
A admixtures in the nuclear wave functions [39], are neglected. However, in 
Ref. [27] it was shown that these approximations do not influence in any 
significant way the theoretical predictions for the pp weak capture cross section 
once the two-body current model is constrained to fit the GT matrix element 
of tritium. 

The results for A^(0) including axial two-body currents are listed for the 
AV18 and CD-Bonn interactions in Table 4. Predictions for this quantity with 
other modern interactions are expected to be similar. Thus, the model depen- 
dence and theoretical uncertainty appear to be at the level of a few parts in a 
thousand, much smaller than the estimate given in Ref. [23]. 

More recently [45], the pp (and hep, see below) astrophysical rates have 
been calculated using wave functions obtained from solutions of the Schrodinger 
equation with realistic interactions, and axial two-body currents derived from 
chiral effective field theory (EFT), as discussed at the end of Sec. 4.2. In 
Ref. [45], the pp 5- factor is the physical observable chosen to test the en- 
ergy/momentum cutoff dependence. It has been demonstrated that this ob- 
servable is independent of the cutoff as long as its value is within a physically 
reasonable range. In particular, this work gives for the ratio of the two-body 
to one-body contributions the value (0.86 ± 0.05)%, in line with the results of 
the conventional theory reported above. 

6.3 The d{n, and d{p, ^^He radiative captures 

The theoretical description of the A=3 and 4 radiative and weak cap- 
ture processes constitutes a challenging problem for nuclear few-body the- 
ory [35,33,39,87,88]. The difficulty comes about because the ^4=3 and 4 bound 
states are approximate eigenstates of the Mi and GT one-body operators [89]. 
This property would be exact if these states were to consist only of a sym- 
metric S-wave component, which accounts for over 90% (80%) of the A=3 
(4) bound states. It is spoiled by the presence of D-state components induced 
by tensor interactions. Consequently, the one-body matrix elements of these 
(largely) Mi- or GT-induced processes are suppressed due to orthogonahty 
between the initial and final states. As a result, the associated cross sections 
become extremely sensitive to contributions from small components in the 
wave functions and two-body electroweak current operators. 

The extensive experimental data for the d{n,'j)^B. and (i(p, 7)^He reac- 
tions include both total cross sections and spin polarization observables at 
several center-of-mass energies. The most complete comparison between the- 
ory and experiment for the (i(n,7)"^H and d{p,'~f)^}ie radiative captures below 
the deuteron breakup threshold was performed using CHH wave functions ob- 
tained from the AV18/UIX interaction [40,76], with a nuclear electromagnetic 
current operator that includes both one- and two-body contributions. The 
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Fig. 2. The calculated S'-factor [76] for d{p,-ffRe obtained using AV18/UIX CHH 
wave functions and one-body only (dashed line) or both one- and two-body (solid 
line) currents. Experimental data are taken from Refs. [94,95,96,97]. 

proton radiative capture calculation has been extended above the deuteron 
breakup threshold [90,91,92], but we consider only the low-energy regime. 

The theoretical prediction of Ref. [40] for the total cross section cr-r of ther- 
mal d{n, 7)^H is 0.578 for the full and 0.229 mb for the impulse approximation 
(lA) current operators. The experimental value is ctt = 0.508 ±0.015 mb [93]. 
Therefore, the "full" result, while clearly an improvement over the lA, still 
exceeds the experimental value by 14%. The origin of this discrepancy is puz- 
zling, particularly in view of the fact that the astrophysical S'-factor for the 
related dp radiative capture at zero energy is calculated to be within 1% of 
the measured value (see below). 

The calculated S'-factor [76] in the energy range 0-50 keV is compared 
with the experimental data of Refs. [94,95,96,97] in Fig. 2. The solid curve 
represents the "full" result, while the dashed line is the result in lA. The 
agreement between theoretical predictions and experimental data, especially 
the very recent LUNA data [94], is excellent. The two-body contributions 
to the nuclear electromagnetic current operator play a very important role 
in this agreement. The "full" (lA) theoretical value for S'(O) has been found 
to be 0.219 (0.162) eV-b, in excellent agreement with the LUNA result of 
0.216 ± 0.010 eV-b. The theoretical and experimental S'-factors also agree 
nicely at higher energies, as shown in Refs. [90,92] in the energy range 0-2 
and 0-3.33 MeV. 

6.4 The hen and hep reactions 

The total cross section of the ^He(n, 7)^He (hen) reaction is experimentally 
known [98], while the astrophysical S'-factor of the ^He(p, e"*" z/e)^He (hep) reac- 
tion cannot be directly measured in the energy range relevant for solar fusion. 
The most recent ab initio calculation of the hen reaction was performed more 
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than a decade ago in Refs. [39,87], using VMC wave functions obtained with 
the older AV14/UVIII Hamiltonian model. The total cross section at thermal 
neutron energy was found to be 85.9 /xb when one- and two-body contribu- 
tions to the electromagnetic current are included, and 5.65 //b in lA. The most 
recent experimental determination gives 55±3 ^b [98]. The discrepancy be- 
tween theory and experiment could indicate that: i) the VMC wave functions 
obtained with the AV14/UVIII interaction model arc not sufficiently accurate; 
and/or ii) the model for the electromagnetic current needs to be refined. Work 
currently under way to construct more accurate CHH wave functions for the 
four-nucleon states involved in the hen reaction will help to clarify this issue. 

The Super-Kamiokande (SK) measurements of the energy spectrum of 
electrons recoiling from scattering with solar neutrinos [99,100,101,102] have 
brought new interest in the hep reaction. The first SK results for this spectrum 
showed an apparent excess of events in the highest-energy bin [100] that could 
be explained by increasing by a factor of ~ 17 the Standard Solar Model 
(SSM98) prediction [103] of the hep neutrino flux. The most recent results 
from SK are essentially the same after taking the large mixing angle (LMA) 
neutrino mixing parameters into account. However, when the quasi- vacuum 
solution is used, the inferred flux at the earth is reduced by a factor of 4. The 
hep neutrino flux is directly proportional to the 5"- factor at zero energy S{0), 
and the SSM98 estimate for the flux was based on the calculation of Ref. [39]. 
This calculation used AV14/UVIII VMC wave functions, retained only the 
p^Hc ^Si channel, neglected the dependence on the momentum transfer of 
the lepton pair, and used poorly-constrained two-body contributions to the 
nuclear axial current operator. 

Because of the signiflcant progress made with the CHH method in the 
description of the four-nucleon bound and continuum states [47,55], the hep 
reaction has been re-examined in Refs. [35,88] and Ref. [45]. All S- and P- 
wave capture states are calculated, using AV18/UIX interaction and retaining 
the g-dependence of the nuclear weak transition operator. In Refs. [35,88] the 
TCO model for the weak current is used, while in Ref. [45] the EFT model 
is used, as discussed in Sec. 4.2. The momentum cutoff A was varied between 
500 and 800 MeV. The results for iS'(O) summarized in Table 5 show that 
the P-wave capture channels are very important and give about 33% of the 
calculated ^'-factor. Contributions from the D-wave channels are expected to 
be small [35,88]. The A dependence in the ^5*1 channel is the result of the large 
cancellation between the one-body and the two-body contributions found in 
Refs. [35,39,88]; such dependence is much smaller in all other channels. 

The energy dependence of the hep ^'-factor was found in Refs. [35,88] to 
be rather weak at 0, 5 and 10 keV. In fact, the calculated ^'-factor at 10 keV, 
close to the hep solar Gamow peak, is 10.1x10"^° keV-b, only 4% larger than 
that at keV. In Ref. [35], iS'(O) was also calculated using the AV18 and the 
AV14/UVIII interaction models. The AV14/UVIII (AV18) resufi is 4% (20%) 
larger than that for AV18/UIX, demonstrating the need to use a Hamiltonian 
that accurately reproduces the properties of the three- and four-nucleon bound 
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Table 5 

Contributions to ^(O) for hep from individual partial waves in units of 10"^'^ keV-b. 





A (MeV) 


'So 


'Si 


'Po 


'Pi 


'Pi 


'P2 


Total 




500 


0.02 


7.00 


0.67 


0.85 


0.34 


1.06 


9.95 


Ref. [45] 


600 


0.02 


6.37 


0.66 


0.88 


0.34 


1.06 


9.37 




800 


0.02 


4.30 


0.66 


0.91 


0.34 


1.06 


7.32 


Ref. [35,88] 




0.02 


6.38 


0.82 


1.00 


0.30 


0.97 


9.64 



and scattering states. 

Following the results of Refs. [35,88], the Standard Solar Model predic- 
tion for the hep neutrino flux has been recently readjusted to be 9.3 xlO^ 
cm^^s"^ [104], 4.4 times larger than the SSM98 prediction, but stiU about 4 
times smaller than the latest SK result based on LMA [102]. 

6.5 The ci(Q;,7)^Li radiative capture 

The reaction d{a, j)^Li is isoscalar, so direct Ei transitions are suppressed. 
The Ml transition from S-wave initial states is also suppressed because the 
ground state is an approximate eigenstate of the Mi operator. As a result, 
the reaction at most energies between 0.2 and 3 MeV is almost purely E2 
in nature, arising from the D waves. The cross section is in the nanobarn 
range and difficult to measure. Nonetheless, there are two sets of cross section 
measurements [105,106], and one set of capture cross sections inferred from 
the breakup process ^"^^Pb -|-^Li — > ^°^Pb + a + d [46]. This small cross section 
helps preclude significant production of ^Li in the big bang. 

Wc show our computed cross sections [64] with the data in Fig. 3. These 
calculations include all initial states with L < 3 and the Mi, M2, Ei, E2, and 
E3 operators. The two sets of theory curves correspond to two different choices 
of V^sj to generate the inter-cluster correlation ip^d'^ (''^otd) of Eq.(16). There 
is very little dependence on Vlsj, probably because it is well-constrained in 
the D waves by low-lying resonances. 

The total cross sections provide a good match to the data at the narrow 
3+ resonance at 711 keV, partly because Vlsj was adjusted to produce the 
location and width of this resonance accurately [107] in the scattering phase 
shifts. At energies between the resonance and 3 MeV, the energy dependence 
of the cross section provides a good fit to the Robertson data [105]. Beyond 
3 MeV, the model fails to reproduce the data, probably because of omitted 
couplings to channels of nonzero isospin. To match the normalization of the 
data, our curve would have to be renormalized by a factor 0.85, and the 
residuals of the data with respect to the renormalized curve are of the order of 
10% even for the higher- quality of the two data sets reported in Ref. [105]. We 
conclude that the theoretical model has come out quite well. This probably 
reflects a good asymptotic normalization coefficient for projection of the ^Li 
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Fig. 3. S'-factor for the process d{a, 7)^Li, showing the variation in theoretical curves 
due to choice of Vlsj- See text for a discussion of the Robertson Ei [105] points. 



ground state onto ad cluster states. 

In the region below 500 keV, there are only indirect cross sections (and a 
rather high upper limit [108]) for comparison, inferred from the process ^°'^Pb + 
^Li — i> ^°^Pb + a + d [109]. This energy range is very interesting, because at 
about 250 keV, we compute the isoscalar Ei transition to be half of the total 
cross section. It becomes more important at lower energies, as the centrifugal 
barrier disfavors the E2 contributions arising from the D waves. We do not 
find the same energy dependence as the indirect measurements, despite our 
inclusion of terms up to third order in the long-wavelength approximation 
(LWA) El operator that are important because of the isoscalar transition. The 
energy dependence of the data suggests S-wave capture, but we compute the 
cross section for E2 capture from the S wave to be two orders of magnitude 
smaller. The angular distributions of the breakup data are consistent with 
purely E2 transitions, but the technique also gives much less weight to Ei 
transitions than to E2. 

Although the Ei operator makes a small contribution to the total cross 
section above 1 MeV, it interferes with the dominant E2 amplitudes to produce 
a forward-backward asymmetry in the distribution of emitted photons. It is 
a persistent difficulty of capture models that they agree with the laboratory 
data on neither the sign nor the magnitude of the asymmetry[105,110]. Ours 
also disagrees with the data, despite the inclusion of several corrections to 
the LWA El operator, because the leading effect remains the center-of-energy 
effect arising from the difference in mass per nucleon between the two initial- 
state nuclei. We note, though, that the upper limits on asymmetry at 1.0 
and 1.63 MeV are much smaller than one would expect from extrapolating 
to lower energies the Ei contribution measured at 2 MeV in essentially any 
capture model. This suggests that the shortcoming may be in the data. 
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Fig. 4. S-factors for a captures forming ^Li and ^Be, showing the spread in the 
computed values due to choice of Vlsj] values for ^H(a,7)^Li are multiplied by 2. 

6.6 The ^He(Q;, 7) ''Be and '^H(a, 7)'^Li radiative captures 

The a-capture processes ^He(a,7)''Be and ^H(a,7)^Li are closely related, 
and have also been computed by the VMC method [65]. We used the same 
phenomenological potentials Vlsj of Sec. 5.4 to describe the continuum in 
both systems, except that the charges in the Coulomb terms and the reduced 
masses were adjusted to match each system. Both processes are Ei transitions, 
originating from S waves at the lowest energies and also partly from D waves 
at the higher end of the energy range shown in Fig. 4. We computed operators 
for the same L range and multipolarities as for d{a, 'j)^Li and found that 
no transitions besides Ei matter above the 0.1% level for either reaction. 
Corrections to the LWA Ei operator are also negligible (though the center- 
of-energy correction, omitted because it is expensive to compute, should be 
on the order of 3% in these cross sections). For both reactions, capture may 
occur into either of two bound states, with a branching ratio of roughly 0.4. 

The computed cross sections [65] for the process ^H(a, 7)^Li are shown and 
compared with the data [111,112,113,114] in Fig. 4. By far the most precise 
experiment to measure this cross section is that of Brune et al. [115] (upward- 
pointing open triangles in Fig. 4). Our results for this reaction have a disper- 
sion of about ±5% among calculations for different choices of cluster-cluster 
potential Vlsj- This dispersion, wider than for (i(a,7)^Li or ^He(a, 7) ''Be, 
probably reflects a substantial contribution to the capture from small a^H sep- 
arations. Taking account of this dispersion and of a roughly 10% normalization 
uncertainty from our Monte Carlo integration, we are in excellent agreement 
with all the modern data. The ratio of the cross section for capture into the 
excited state of ^Li* to that for the ground state is almost independent of Vlsj- 
While the calculation reproduces the approximate absence of energy depen- 
dence seen in the branching ratio data, the predicted ratio itself is about 15% 
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below the data. (The branching data have a 4% normahzation uncertainty.) 

The computed energy dependence of the '^H(a,7)''Li total cross section 
agrees reasonably well with the Brune data, for a reduced of 2.5. This 
compares favorably with other calculations in the literature. At the highest- 
energy data points, where the S'-factor is increasing due to a growing contri- 
bution from D-wave capture, our calculation has a slightly shallower energy 
dependence than the data. We also find a significantly shallower slope for the 
S'-factor at zero energy than all but one other theoretical calculation [11]. 

Fig. 4 also shows the computed [65] and measured [116,117,118,119,120,121] 
yS-factors for the process ^He(Q;, 7)^Be. Because of the larger Coulomb barrier, 
this reaction is more peripheral than ^H(q;, 7)^Li, and this is probably why its 
cross section depends less on Vlsj- This is probably also why there is excellent 
agreement in the energy dependence of the cross section both with the labora- 
tory data and with older theoretical models. The cross section normalizations 
of the various experiments are not completely consistent with each other, but 
they all show higher normalizations than do the results of our calculation, 
ranging from a 10% to a 100% discrepancy relative to our curve. While a 10% 
discrepancy is within the range allowed by the combined uncertainties of the 
experiments, of Vlsj, and of our Monte Carlo integration, the more discrepant 
data sets are in clear disagreement with our result. Since the reaction is pe- 
ripheral, and the branching ratio for capture into the two bound states is well 
reproduced, this suggests that the wave functions for both ^Be bound states 
have incorrect values for their asymptotic normalizations in the L = 1 a ^He 
channels, and by about the same factor. Because of the discrepancies among 
experimental data, this reaction is a case where a microscopic model that 
can reliably predict the cross section normalization has the potential to shed 
light on the nature of the experimental difficulties and to be very useful to 
astrophysics. Future calculations with wave functions that are more accurate 
solutions for the underlying two- and three-nucleon interactions (variational 
or GFMC continuum states and GFMC bound states) should achieve that 
goal. 

6.7 ^ Be weak decay 

The ^Be nucleus decays by electron (e-)capture to the ground state of ^Li 
and to its first-excited state at 0.478 MeV. These decays show up in the solar 
neutrino spectrum as two sharp lines at 0.862 and 0.384 MeV, respectively. The 
two ground states are both spin-| states, while the ^Li* excited state is spin-|. 
Consequently there are both Fermi and Gamow- Teller (GT) matrix elements 
to ^Li, but only GT matrix elements to ^Li*. These have been calculated using 
both Type 1 and Type 11 VMC wave functions as described in Sec. 5.2, in 
impulse approximation (lA) and with one-body relativistic, two-body meson- 
exchange-current, and isobar contributions [36]. 

The Fermi matrix element is F = -J'iJf + 1 for Type I wave functions 
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Table 6 



Gamov- Teller matrix elements for A=3,7 nuclei. 


CHH 


VMC(I) 


VMC(II) 




lA Full 


lA Full 


lA Full 


Expt 


^^He 1.597 1.658 


1.602 




1.658 


^Be ^^Li 


2.345(3) 2.419(5) 


2.367(2) 2.433(5) 


2.599 


■^Be ^^Li* 


2.142(2) 2.200(3) 


2.141(1) 2.205(3) 


2.323 



due to isospin symmetry, but is slightly less, F = —1.999, for Type II wave 
functions which include long-range Coulomb correlations. The results for GT 
matrix elements are shown in Table 6, where we also show the matrix ele- 
ments used to fix terms beyond lA. (The A — 7 lA results have been updated 
from Ref.[36] to reflect recent improvements in the VMC wave functions.) The 
lA matrix elements have been calculated with both CHH and VMC wave 
functions with agreement to better than 0.5%. The relativistic, mesonic, and 
isobaric corrections total about 4% in ^H, and 3% in ^Be. The two types of 
VMC wave functions for A = 7 give the same results within 1%. Compared to 
experiment, however, the total GT matrix elements are too small by 7% for 
the ground state and 5% for the excited state; the corresponding half-life is 
calculated to be 62.3 days in lA, and 58.2 days in total, compared to the exper- 
imental 53.22±0.06 days. The branching ratio to these two states is calculated 
to be 10.20% in lA, and 10.33% in total, compared to an experimental value 
of 10.44±.04%. These calculations should be repeated with the more precise 
GFMC wave functions, and with an improved Hamiltonian such as AV18/IL2 
that does not under bind the A — 7 nuclei. 

7 Conclusion and prospects for the future 

Recent advances in the descriptions of nuclear Hamiltonians, wave func- 
tions, and electroweak currents - as well as computer speed - may be fruitfully 
applied to several reaction and decay processes of interest for astrophysics. Cal- 
culations based on realistic NN interactions and currents may be relied upon 
for a more fundamental understanding of the nuclear systems than previous 
treatments. In the case of the weak captures, experimental fitting or verifica- 
tion of cross sections is not possible, so it is absolutely crucial to have the most 
fundamental theoretical understanding possible. The radiative captures arc in 
many cases observable in the laboratory, but calculations are still necessary 
as tests of wave functions and currents, as guides for extrapolating laboratory 
data to lower energies, and (eventually) as independent sources of information 
on the absolute cross section. 

There remains considerable opportunity for future efforts in this area. For 
captures producing p-shell nuclei, work is under way to develop the VMC 
technique for continuum states so that they can be derived directly from the 
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underlying NN interaction. Future improvements will include the use of the 
new Illinois 3iV interactions to generate wave functions, as well as the devel- 
opment of methods to compute matrix elements from the (essentially exact) 
GFMC wave functions. More direct extraction of the asymptotic normaliza- 
tion coefficients from the quantum Monte Carlo wave functions may in the 
future prove to be very useful for direct captures. 

More reactions will also be accessible in the future, partly because of im- 
provements in algorithms and in computer speed. Ongoing work will extend 
theoretical calculations of the reaction d{p, 7)^He to higher energy. The reac- 
tion ^Be{p, 7)^B is an important problem in light of new laboratory data and 
the need to extrapolate them to solar energies. The isospin mirror process, 
'^Li(r2, 7)^Li has been extensively measured [122] and will provide an interest- 
ing test of the method. Also in progress is a study of the weak decay of ®B, 
which is much like the hep process in the presence of a spectator a-particle. 
Some longer-term prospects are ^B(p, 7)^0 and ^He(Q;n, 7)^Be reactions that 
may provide bridges across the A = 5 and A = 8 stability gaps in some as- 
trophysical environments. The latter will most likely have to be treated as a 
neutron capture on the lowest 0"*" state of *Be. 

As more is learned about solar neutrinos, the big bang, and other areas 
of nuclear astrophysics, the demand for greater precision in nuclear inputs 
to astrophysical calculations is increasing. Ongoing advances in theoretical 
treatments of reactions and decays among the light nuclei will continue to 
help meet this demand. 
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